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Abstract 

0\ . 

Recent measurements of B° — > tttt decays reveal two features that are in conflict with 
conventional calculations: the channel B°(B°) — > n + ir~ shows a large direct CP- violating 
asymmetry, and the channel B°(B°) — > 7t°tt has an unexpectedly high branching ratio. 
\ We show that both features can be understood in terms of strong-interaction mixing of 

f^S ■ tttt and DD channels in the isospin-zero state, an effect that is important because of the 

large experimentally observed ratio T(B°/W -> D+ D-)/T(B° fW -> tt+tt") w 50. Our 
dynamical model correlates the branching ratios and the CP-violating parameters C and S, 
for the decays B°{W) -> tt+tt", B°(W) -► tt tt°, B°(B°) -► £>+£>- and ^"(BO) -► 

The Belle collaboration has presented new data which support their original evidence 
[2] for large direct CP violation in the decays B°(B°) — > tt + tt~, the asymmetry parameter C 
(= —A) being measured to be C = —0.58 ± 0.15 ± 0.07. In a related development, both the 
Babar [3] and Belle [I] collaborations have reported a sizable branching ratio for the decay 
B°(W) -» 7T°7T°, with an average value Br(B°/W 7r 7r°) = (1.9 ± 0.6) x 10~ 6 . Both of these 
observations are unexpectedly large from the standpoint of conventional calculations El Ej 
based on a short-distance, effective weak Hamiltonian and the assumption of factorization of 
products of currents in matrix elements for physical hadron states. In this paper, we carry out 
a calculation based upon the idea (Hj of final-state interactions involving the mixing of 7T7T and 
DD channels. This dynamics provides a natural, correlated explanation of the new experimental 
facts, and leads to several further predictions. 

To fix notation, we write the three B — > tttt amplitudes as 

A(W^tt + tt-) = N(X uai + X c a p ) 

A(B~v^tt tt ) = N(X u a 2 -X c a p )/V2 (1) 
A(B~ -» tt-tt ) = NX u ( ai + a 2 )/V2 



Here a±, a 2 , a p are, in general, complex numbers and N is a positive normalization factor. The 
parameters X u and A c are CKM factors, defined as X u = V u bV* d , X c = V&V^, with magnitudes 

* E-mail address: sehgal@physik.rwth-aachen.de (L. M. Sehgal) 



1 



|A U | = 3.6 x 10~ 3 , |A C | S 8.8 x 10" 3 and phases given by X u = IX^e'^, X c = -\X C \, with 7 ra 60° 
[S]. The amplitudes in Eq. Q are defined so that their absolute square gives the branching 
ratio, and they satisfy the isospin relation ^U] 



T2 A{B " 



Tr+vr-) + A(B° -» vrV) = A(B~ -» 7r~7r' 



(2) 



From the results of the models discussed in El El) the parameters appearing in Eq. Q 
have the following rough representation. The constants a±, a>2, a p are approximately real (to 
within a few degrees), with magnitudes a\ ~ 1.0, 02 ~ 0.2, a p ~ —0.1. The normalization factor 
is N ~ 0.75; it is here fixed by the empirical branching ratio for B~ 



TT 7T°. 



The fact that 

the parameters a\, 0,2, a p are nearly real implies immediately that there is very little direct 
CP-violating asymmetry between B° —* tt + tt~ and B° — ► 7r + 7r~, as well as in the channels 7r°7r° 
and 7r ± 7r°. Furthermore, the absolute branching ratios following from the above parametrization 
are as follows (with experimental values given in parentheses): 



Br( J B ± - 
Br(B°/W - 
Bt{B°/W 



7T 7T 



TT TT ) 
-> 7T 7T J 



5.3 x 10- 
9.2 x 10" 
0.2 x 10" 



[exp. (5.3 ± 0.8) x 10" 6 ] 
[exp. (4.6 ± 0.4) x 10~ 6 ] 
[exp. (1.9 ± 0.6) x 10~ 6 ] 



(3) 



The most striking feature is the strong enhancement of the rate compared to this model 

expectation. 

It was pointed out in Ref. jS] that the CP- violating asymmetries and branching ratios in 
the B — * tttx system would be strongly affected by final-state interactions involving the mixing 
of the 7T7T and DD channels in the isospin 1 = state, as a consequence of the large ratio of 



TT 



11 



cb 



w, 



ub 



26 expected in the 



partial decay widths r(P° -» D + D-)/T(B° -► vr H 
Bauer-Stech-Wirbel model 5 . A large ratio has now been confirmed by the Belle measurement 
jllj of the branching ratio Br(P°/P° — ► D + D~) = 2.5 x 10~ 4 , which is about 50 times larger 
than Br(P°/P° — ► tt + tt~). This fact gives new urgency to an investigation of tttt «-> DD mixing 
as a way of resolving the puzzling observations in B — > tttt decays. 

The tttt system exists in the states I = or / = 2, while the DD system has / = or / = 1. 
Mixing can occur between the isospin-zero states 



F7T/0 
1^)0 



2 | + ~ 

3 1 

r 



\D+D~) + \D°D°) 



1|7T 
0" 



(4) 



By contrast, the 1 = 2 tttt state and the 1 = 1 DD state, given by 



F 7r /2 



\DD)i 



7T + TT ) 



\D + D~ 



TT TT ) 



\D°D°) 



(5) 
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are unaffected by mixing. The physical decay amplitudes of B° to the above four states are 



A (0) 



AO) 

DD 

Ai) 

DD 



D- 



A 



D+D- 



3 

+ A 
-A 



(6) 



D°D° 



D°D° 



These physical decay amplitudes are related to the "bare" amplitudes calculated in the absence 
of final-state interactions, i.e. with no mixing, which we denote by A: 



(7) 



Here S denotes the strong-interaction S matrix connecting the isospin-zero states \tt7t)o and 
\DD)q which can be written generally as 1 



4 (0) \ 

AO) 

DD ' 


= s 1 ^ 


AO) \ 

AO) 

DD > 


J '"KIT 


- AW 




DD 


~ DD 





s 



cos26e i2Sl ism26e i( - Sl+5 ^ 
ism26e i( - Sl+s ^ cos26e i252 



(8) 



where 9 is a mixing angle, and 5\ and 62 are the strong-interaction phase shifts for the elastic 
scattering of -kit and DD systems in the 1 = state, at = Mb- For any choice of these 
three parameters, the matrix can be calculated numerically, and the set of four equations 



and A DO jyj i n terms of the 



(J7J) solved to obtain the physical amplitudes A n + n -, A n o n o, A D + D - 
bare amplitudes. The bare amplitudes are identified with those calculated in the factorization 
model El El j which we list below 



A j 

■™"ir + n~ 

A-d+D- 
^D°D° 



iV(A u ai + X c a p ) 
N(X u a 2 - X c a p )/V2 
N'Xcax 




(9) 



where the first two equations are as in Eq. (Q) , and the factor N' is determined from the empirical 
P] branching ratio Br(B°/B°~ -> D+D~) = iV /2 |A c | 2 a 2 = 2.5 x 10 -4 to be N' = 1.79. 

In order to show, in a transparent way, how the mixing mechanism gives rise to large direct 



7r°7r°, we 



CP violation in B° — > ir + ir , as well as an enhanced branching ratio for B° 

1 The two-channel S'-matrix has been discussed, in particular in |12II13| . The prescription is given in pjll!2|. 
An alternative prescription, using i [1 + S] in place of has been discussed by Kamal |14|. and was used in 
Ref. 0. 
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consider, for illustration, the case where the elastic phases 8± and 82 in the S matrix (Eq. (jSJ)) 
are neglected, so that £2 may be written as 



cos 9 i sin 9 
i sin 9 cos 9 



(10) 



The amplitudes A n + n - and A n o n o for B° decay are then given by 

1 + 2 cos 



- A^ 



r- 1 — COS 9 ~ 

+ V2 A^a^o 



1 



1^-= [A D +D- +^ D o^o) 



A^o^o 



V2 l - COS °A^ 



+ 



2 + cos 9 



A n o n o 



-jsini 



D+D- 



+ A 



D°D°) 



(11) 



Clearly for 9 = 0, the physical amplitudes reduce to the bare amplitudes. Inserting the bare 
amplitudes from Eq. @, we can rewrite A 1T + 1T - and A^o^o as linear combinations of \ u and A c : 



A r , 



N 

N 
72 



.1 + 2 cos 9 1 — cos i 1 
a u <J a\ H 02 \ + A c (a p cos # + a ri 



x . 2(1- cos 0) 
A« 4 — 7, -a-i + 



3 

2 + cos # 



:i2) 



02 ^ — A c (a p cos + a m ) 



where 



.1 . ,N> 
z— = sm#— — ai 



(13) 



Note that the isospin relation in Eq. © continues to be fulfilled. The important new feature 
of the amplitudes in Eq. (|12|) is the appearance of the imaginary term a m in the coefficent of 
A c , in striking contrast to the real term a p . The imaginary nature of this dynamical term is an 
inescapable consequence of 5-matrix unitarity, which enforces the factor i in the off-diagonal 
matrix element in Eq. (|1()|) . The term a m , given in Eq. (|13|) . has a magnitude \a m \ ~ 1.39 sin 9, 
and dominates the term a p cos 9 even for a modest mixing angle ~ 0.1. We will now show that the 
mixing term a m has profound consequences for direct CP violation in the decays B° — > 7r + 7r~, 
and for the branching ratio of the channel B° —* 7r°7r°. 



1 C and S Parameters for B 







7T + 7r and B 







The C and 5 parameters derived from the time-dependent asymmetry between B° and B° decays 
into 7r + 7r~ are defined as 



s 4 



1-IA, 



1 + IA+-I 2 
2ImA+_ 
1 + |A + _| 2 



(14) 
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where 



with 



q 
p 



q A(B° - 
pA(B° - 

e~ i2f} , 2/3 



TT ' TT J 
TT + TT~) 

3 45°. 



C-| is the parameter for direct CP violation (i.e. 

Using the amplitude A n + n - in Eq. (|T2*|) . we obtain 

A 



A, 



-i2/3 



1+2 cos I 1— cos 

o Ul -|- 



7r+7r-)/A(B° 
;.» 02) + A c (a p cos + a m ) 



TT ' TT 



A? 



1+2 COS 9 1 1; 

3 ai ^~ 



«2 ) + A* (a p cos 9 + a m ) 



(15) 
(16) 

(17) 



The asymmetry parameters C_| and calculated from the above expression are plotted as 

functions of 9 in Fig. ^ Good approximate agreement with data is obtained for 9 ~ 0.2 (see 
Tabled where we also list Coo an d Sqq). We note that the amplitudes in Eq. ()12|) have been 
derived from the matrix S2 in Eq. (|1U|) . which was obtained from (JHJ) by neglecting the phase 
shifts Si and 62- We have also explored numerically 5 matrices with non-zero phases, and 
indicate in Figs. ^ an d[21two examples, obtained with the values <5i = ±10°, 61 + 62 = —30°. 
Table^gives numerical values for a few choices of parameters. In all cases, there is a large direct 
CP violation. 

Discussions of the direct CP-violating parameter C_| are often based on an amplitude for 

B° — > it + 7t~ written in the form 

P ' 

1 7T7T 



A Tl 



-17 



+ 



(18) 



The parametrization in Eq. based on the models El U\, gives \P n7T /T nn \ = 0.24, and 
arg(P 7Tn /T nn ) = 0. The small phase of the "penguin-to-tree" ratio P^/T^ is a generic feature 

of these models, and is responsible for the prediction C_j ~ 0, which is now contradicted by 

data p. In our approach, the role of P^/T^ is played by the ratio 



U P /T n 



I A c I (dp cos 9 + a„ 



IA 



.i( 



1+2 cos 9 1 1— cos ( 
3 Ul 3 



a 2 



(19) 



For a typical value 9 = 0.2, this ratio has the modulus | u p/T" \ ~ 0.77, and a phase arg( "P/T") ~ 
—70°. The difference is a consequence of the term a m in Eq. (|19|). which reflects the physical 
final-state interaction of the tttt system, as implemented in our model through tttt <-> DD mixing. 



2 Branching Ratio for B —> tttt and B — > 7t + tt 

The branching ratios (averaged over B° and B°) may be calculated in our model by taking the 
absolute square of the B° decay amplitudes in Eq. (|12[). and the corresponding amplitudes for 
B° decay. The results are shown in Fig. [21 It is remarkable that the empirical branching ratio 
for B°/B° — > 7r°7r° is accurately reproduced, using the same value 9 ~ 0.2 which accounts for 

the asymmetry parameter We also note that the branching ratio B°/B° — > it + tt~ remains 

close to its bare value, and can be lowered slightly with the introduction of phases 6± and 62- 
Numerical results for Bv(B°/W -» vr vr ) and Br(P°/P° -> vr+vr-) are listed in Tabled 
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3 Branching Ratio for B° -> D°D° 

Since our model treats the tttt and DD states with / = as a coupled system, it also produces 
predictions for branching ratios and asymmetry parameters in B° — > D + D~ and B° —* 
The amplitudes after mixing are 



A 



D+D- 



isin6»y^ (y/2A n+n - - A^o 



+{cos6 + l)A D+D - + (cos9-l)A D0W (20) 
1 " 



I sin i 



+(cos - l)i D+D - + (cos + l)A D0W 



Of particular interest is the branching ratio for B°/B° — > D°D°, since it vanishes at the level 
of the bare amplitude (A^-g^ = 0), and is induced by mixing with the 7T7t system. For 9 = 0.2, 
ignoring the phases Si, $2, our model predicts 

Bt(B°/W -> Z>°!d°) = 1.45 x 10~ 7 (21) 

(At this low level, one must assume that other sources of final-state interaction or a non-zero 
bare amplitude could raise this branching ratio further.) Direct CP violation follows from A D0 -^ 
in Eq. (J^UJ): C D0 ^ = -0.50 for 9 = 0.2. Direct CP violation in D + D~ (and in vr^vr ) is small, 
because these decays are dominated by a single amplitude. There is little mixing in A^+d- in 
Eq. (|20|) (and none in the 1 = 2 amplitude for 7r _ 7r°). 

To conclude, we have demonstrated a mechanism of final-state interactions among physical 

hadrons in B° — ► tttt decays which predicts a large direct CP-violating parameter C_| The 

same mechanism enhances the theoretical prediction for the branching ratio of B° /B° — > 
to the experimentally observed level. Predictions are made for the C and S parameters of 
B°(B°) — > 7T°7r decays, and for the branching ratio of B°/B° — > D°D°. The model makes 
essential use of the large empirical ratio r(£°/#° -> L>+ £T) /T(B° fW -> vr+vr") fs 50. Its 
success in the present context leads to the expectation that sizable direct CP violation could be 
observed in other charmless B decays, in which an amplitude of order \ u receives a dynamical 
contribution proportional to A c , through mixing with a channel possessing a large branching 
ratio. The resulting amplitude contains two pieces which are comparable in magnitude and 
have different weak-interaction and strong-interaction phases. We have treated earlier |15] the 
charged-particle decays B^ — > rjir^ (and B^ — * rj'ir^), which are influenced by mixing with the 
channel B^ — * r]^^ 1 , and have predicted significant direct CP violation. Evidence for a sizable 
violation in B^ — * rjtr^ has indeed been reported in one experiment the first ever seen in 
a charged-particle decay. 
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Observable 


No 
mixing 


With mixing 


Data 


9 = 0.2 

51 = 0° 

5 2 = 0° 


e = 0.17 

51 = -10° 

5 2 = -20° 


9 = 0.2 

St = 10° 
5 2 = -40° 


C+- 
S+- 

Br(B°/W -> ^°vr ) 
Bt(B°/W ^vr+vr-) 


±0.00 
-0.60 

0.2 

9.3 


-0.65 
-0.63 
1.8 
12.2 


-0.66 
-0.55 
1.7 
10.5 


-0.81 
-0.40 

1.6 

9.9 


-0.58 ± 0.15 ± 0.07 (Belle 0) 
-0.30 ± 0.25 ± 0.04 (Babar [IZI) 
-1.00 ± 0.21 ± 0.07 (Belle 0) 
±0.02 ± 0.34 ± 0.05 (Babar [TZ]) 
1.7 ±0.6 ±0.2 (Belle 0) 
2.1 ±0.6 ±0.3 (Babar 0) 
4.4 ± 0.6 ± 0.3 (Belle HBJ) 
4.7 ±0.6 ±0.2 (Babar [T7|) 
^-\i-ol (CLEO H2|) 


Coo 
<?oo 


±0.00 
±0.73 


±0.48 
-0.65 


±0.51 
-0.78 


±0.56 
-0.49 





Table 1: Observables for different mixing 
units of 10" 6 . 



angles 9 and 



strong-inter; 



action phases Si and S 2 ■ All branching ratios are given in 




9 



6 x KT 6 ] 




Figure 2: Average branching ratios for B°/B° — > 7r°7r° and B°/B° — ► 7r + 7r . Full line is for 
<5i = <5 2 = 0°, dotted line for <5i = -10°, <S 2 = -20°, dashed line for Si = 10°, 5 2 = -40°. 
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